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POSSIBLE  STATISTICS  OF  TWO  COUPLED  RANDOM  FIELDS: 
APPLICATION  TO  PASSIVE  SCALAR 

B.  DUBRULLE’*  AND  GUO  WEI  HE+ 


Abstract.  We  use  the  relativity  postulate  of  scale  invariance  to  derive  the  similarity  transformations 
between  two  coupled  scale-invariant  random  fields  at  different  scales.  We  find  the  equations  leading  to 
the  scaling  exponents.  This  formulation  is  applied  to  the  case  of  passive  scalars  advected  i)  by  a  random 
Gaussian  velocity  field;  and  ii)  by  a  turbulent  velocity  field.  In  the  Gaussian  case,  we  show  that  the  passive 
scalar  increments  follow  a  log-Levy  distribution  generalizing  Kraichnan’s  solution  and,  in  an  appropriate 
limit,  a  log-normal  distribution.  In  the  turbulent  case,  we  show  that  when  the  velocity  increments  follow 
a  log-Poisson  statistics,  the  passive  scalar  increments  follow  a  statistics  close  to  log-Poisson.  This  result 
explains  the  experimental  observations  of  Ruiz  et  al.  about  the  temperature  increments. 

Key  words,  scale  invariance,  scale  transformation,  scaling,  turbulent  scalar  field 

Subject  classification.  Fluid  Mechanics 

1.  Introduction.  Scale  invariance  [1]  usually  refers  to  physical  fields  keeping  the  same  properties  at 
different  scales.  Therefore,  scale  invariance  is  often  associated  with  the  existence  of  power  scaling  laws.  In 
this  case,  the  central  characteristic  quantities  of  the  system  are  the  scaling  exponents.  When  a  random  field 
is  statistically  scale  invariant,  each  of  its  structure  functions  follows  a  power  law,  and  there  is  a  whole  family 
of  exponents  associated  with  its  statistics. 

To  understand  the  underlying  properties  associated  with  the  scaling  exponents  one  can  follow  two  paths: 
first,  try  to  measure  or  compute  the  scaling  exponents  with  as  high  a  precision  as  possible.  This  is  the  path 
usually  followed  in  critical  phenomena,  where  scaling  exponents  are  few,  and  appear  to  be  universal,  due  to 
some  underlying  conformal  symmetry  [2].  In  turbulence,  this  path  seems  more  difficult  to  follow,  because  the 
number  of  scaling  exponents  to  measure  or  compute,  in  order  to  know  the  statistics  of  the  field,  is  infinite. 
Therefore,  parallel  to  the  measurement  or  actual  computation  path,  there  seem  to  be  the  necessity  of  a  more 
formal  path,  devoted  to  understand  link  between  the  various  scaling  exponents.  That  way,  the  measure 
or  computation  of  only  a  few  scaling  exponents  would  hopefully  enable  the  unraveling  of  the  statistical 
nature  of  the  system.  A  natural  way  to  establish  link  between  scaling  exponents  is  to  use  the  symmetry 
responsible  for  the  existence  of  the  scaling  exponents,  namely  the  scale  invariance.  This  was  the  basis  of 
conformal  theory  [2]  or  multifractal  theory  [3].  More  recently,  Dubrulle  et  aJ.  [1,  4,  5,  6]  found  that  new 
constraints  on  scaling  exponents  stemming  from  scale  invariance  could  be  found  via  an  analogy  between 
scale  invariance  and  relativistic  dynamics  inspired  by  a  theory  of  Nottale  [7]  and  of  Pocheau  [8].  This  can 
be  used  to  classify  the  possible  forms  of  scaling  exponents  according  to  the  fixed  points  of  the  similarity 
transformation.  However,  these  results  were  restricted  to  one  random  scale  invariant  field  (ID  case).  In 
many  physical  situations,  the  dynamics  involves  at  least  two  coupled  random  fields,  such  as  transports  of 
passive  scalars  (coupling  of  velocity  and  scalar),  magnetohydrodynamics  (coupling  of  velocity  and  magnetic 
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field)  and  turbulence  (coupling  of  transverse  and  longitudinal  velocity  increments).  The  goal  of  the  present 
paper  is  therefore  to  extend  some  of  these  results  to  the  case  of  two  coupled  scale  invariant  random  fields 
(2D  case)  and  use  the  2D  results  to  explore  scaling  of  passive  scalars. 

The  transport  of  passive  scalars,  for  example  the  temperature  is  governed  by  the  diffusion  equation: 

(1.1)  5t(9  +  (u  •  V)<9  =  k;A(9, 

where  the  solenoidal  velocity  field  u  is  either  taken  to  be  random  or  governed  by  the  Navier-Stokes  equation, 
K  is  molecular  diffusivity.  We  are  interested  in  scaling  of  passive  scalar  structure  functions  for  given  scaling 
of  velocity  structure  functions: 

si{e)  =  mx +£)-  e{x)r)  oc 

(1.2)  S^(e)  =  {{u{x  +  £)-  n{x)r)  oc 

where  ^(n)  and  ^(n)  are  the  scaling  exponents  of  velocity  structure  functions  and  scalar  structure  functions 
respectively.  We  will  use  2D  similarity  transformations  to  derive  these  possible  scaling  exponents. 

This  paper  is  organized  as  follows.  After  useful  preliminaries  regarding  scale  invariance  (Section  2),  the 
2D  similarity  transformation  of  scale  invariant  fields  is  presented  in  Section  3.  The  results  are  applied  to  the 
case  of  the  transport  of  a  passive  scalar  by  a  Gaussian  velocity  field  in  section  4  and  by  a  turbulent  velocity 
field  in  section  5.  We  conclude  the  paper  in  section  6. 

2.  Preliminaries  on  scale  invariance.  Scale  invariance  is  a  commonly  used  concept  but  sometimes, 
rigorous  definitions  are  lacking.  In  this  section,  we  give  a  brief  review  on  the  rigorous  definition  of  scale 
dilation  (or  contraction),  global  and  local  scale  invariance.  Then,  we  discuss  the  existence  of  similarity 
transformation  between  measurements  of  physical  fields  at  different  scales  due  to  scale  invariance. 

Let  ((/?£(x),^^(x))  be  a  two-dimensional  positive  homogeneous  random  field  [9]  and  dependent  on  scale 
We  define  a  family  of  scale  dilation  Sh^h'  (A)  via: 

(2.1)  ^h,h'W  •  ^ 

A  two  dimensional  homogeneous  random  held  is  said  to  be  globally  scale  invariant  or,  simply  scale  invariant, 
if  ((p£(x),‘ip£(x))  and  Sh,h'(X)((p£(x),i/j£(x))  have  the  same  statistical  properties. 

In  log-variables,  T  =  ln£,  X  -  ln((pi),  Y  =  ln(ip£),  the  globally  scale  invariance  amounts  to  a  transla¬ 
tional  invariance.  This  invariance  is  connected  with  the  possibility  to  multiply  any  characteristic  quantity 
by  an  arbitrary  constant,  i.e.  with  the  possibility  to  perform  arbitrary  changes  of  units  in  the  system  [1].  It 
is  a  global  scale  invariance  that  leads  to  the  existence  of  power  laws  of  the  moments  of  random  homogeneous 
fields: 

(2.2) 

Scale  invariance  can  be  also  considered  as  local  [8].  Local  scale  invariance  is  connected  with  the  invariance 
of  the  system  with  respect  to  resolution  changes.  Indeed,  consider  a  discrete  slicing  of  the  scale  space  i  <  L 
in  the  form: 


£i  =  T%,  x  =  r^4, 
ipi  =  A  Vo,  '4>i  =  A  Vo, 


(2.3) 
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where,  T  (or  equivalently  N)  and  A  characterize  the  resolution  of  the  observations,  ipo  and  the  units  of 

the  observations.  By  local  scale  invariance,  we  require  that  the  physical  properties  of  the  system  be  invariant 
under  changes  of  resolution,  i.e.  under  changes  of  N.  In  fact,  changes  of  N  can  be  achieved  by  dilation  of 
F:  r  F",  which  transform  i  into  ia.  It  is  also  a  transformation  which  transforms  li/ij  into  and 

(j)i/(t>j  into  i(l>il(t>j)^^  Therefore: 

A  local  scale  dilation  can  be  seen  as  a  transformation  for  arbitrary  a: 

(2.4)  moT;  i’/^o  i^Phor. 

A  two  dimensional  homogeneous  random  Reid  is  said  to  be  locally  scale  invariant  if  the  physical  Reids 
((p£(x),'ip£(x))  are  statistically  invariant  under  local  scale  dilation. 


Local  scale  invariance  can  be  also  seen  as  a  change  of  rationalized  unit  system  because  it  amounts  to 
multiply  smd  by  suid  Therefore,  global  scale  invariance 

implies  local  scale  invariance. 

When  the  field  {(p£{x),'ipi{x))  is  globally  scale  invariant,  all  these  changes  of  units  can  be  characterized 
by  only  one  number:  the  exponent  of  their  power  law.  This  number  labels  the  system  of  units  and  varies  with 
it.  It  is  not  an  absolute  quantity.  All  these  remarks  have  been  so  far  considered  as  rather  obvious  because 
they  do  not  lead  to  any  specially  useful  information  when  applied  to  deterministic  fields.  For  random  fields, 
however,  <  >  and  <  <j)^  >  can  be  very  different,  and  the  application  of  our  remarks  can  lead  to  very 

interesting  consequences  about  the  exponents.  In  fact,  they  are  not  necessary  additive  anymore  and  can 
obey  more  complicated  composition  law,  which  can  be  determined  only  from  the  requirement  of  local  scale 
invariance,  as  was  shown  in  [6].  We  now  define  the  formalism  which  can  be  used  to  derive  these  composition 
laws  in  the  case  of  two  coupled  fields.  This  is  done  by  a  natural  generalization  of  [1,  6]  which  we  refer  to  for 
more  details. 


3.  Formalism  and  derivation  of  the  composition  laws. 

3.1.  Formalism.  As  noted  previously,  scale  invariance  amounts  to  a  translational  invariance  when 
logarithm  of  scale  and  fields  are  considered.  We  also  would  like  to  introduce  a  field  variable  which  is 
deterministic  but  keeps  track  of  the  possible  different  realizations  of  the  process.  A  good  way  to  reconcile 
these  two  requirements  is  to  introduce  the  log- variables  [6] : 


(3.1) 


{\n{4>il4>o)  {<l>il<hT) 

HMM”) 

_  { Injjji/ipo)  {'ipel'ipoy  ) 


where  K  is  the  basis  of  the  logarithm,  io  is  an  arbitrary  length  unit,  two  arbitrary  field  units  (pos¬ 

sibly  random  and/or  scale  dependent  [1]),  and  n  and  p  are  real  numbers.  We  note  that  In(iF),  p  and  n 
labels  the  system  of  subunits  used  to  measure  the  length  and  fields.  Then,  globally  scale  invariance  implies 
homogeneity  in  these  log-coordinates.  Local  scale  invariance  implies  the  existence  of  a  group  transformation 
connecting  different  measures  of  the  field  at  different  scales.  Formally,  this  means  that  we  can  find  a  simi¬ 
larity  transformation  between  {\ogj^{f/io),Xn,Yp)  and  (log|^/(^/^o)»  Since  £  is  deterministic,  it  is 
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equivalent  to  change  X  or  ^  to  change  T  =  so  we  can  as  well  fix  K  e.g.  K  =  e  and  let  i  vary.  The 

transformation  is  now  between  (T  =  \n{t/io),Xn^Yp)  and  (T'  =  \n{i' /ioj^Xn^^Yp*)-  We  call  it  the  similarity 
transformation. 

The  derivatives  U{n)  —  dXnldT,V(p)  —  dYpjdT  are  analogous  to  the  relative  velocities.  In  our  formal¬ 
ism,  they  are  just  the  local  multifractal  exponents  [6]. 

3.2.  Two-dimensional  similarity  transformation.  The  possible  similarity  transformation  can  be 
derived  provided  the  system  is  both  globally  and  locally  scale  invariant.  Global  scale  invariance  implies 
that  the  transformation  are  linear  [1].  Since  nothing  forbids  a  coupling  between  fields  and  scales,  their  most 
general  representation  is  a  3  x  3  matrix  coupling  X,  Y  and  T.  By  local  scale  invariance,  the  set  of  such 
3x3  matrices  must  form  a  (at  least  semi)  group,  depending  on  two  parameters  U  and  V  [1].  Without  loss 
of  generality,  they  may  be  written  as: 


(3.2) 


A  B  TU 
DEW 
GU  HV  T 


where  A^B^D^E^G^H^U^V  and  T  are  some  unknown  functions  of  U  and  V  to  be  determined.  According 
to  the  transformation  (3.2),  a  field  characterized  by  A  =  Y"  =  0  for  any  T  in  one  reference  field  R  =  (n,p) 
is  transformed  into  X'  =  UT',  Y'  =  VT'  in  another  reference  field  =  {n\p').  This  shows  that  (U,V)  are 
characteristic  numbers  labeling  the  change  from  (n,p)  to  (n',p').  It  is  then  natural  to  adopt  U  and  V  as  the 
group  parameters. 

To  completely  determine  A,  ...,r  as  a  function  of  {U,V),  an  additional  assumption  is  needed.  From 
physical  considerations,  we  choose  to  postulate  that  the  group  of  similarity  transformations  is  commutative. 
This  means  that  it  is  equivalent  to  perform  a  transformation  at  scale  £  then  £'  or  £'  then  £.  This  property  is 
generic  in  the  group  of  scale  dilation  S\{h,  h'),  at  any  given  h  and  h*,  which  motivates  our  choice.  Moreover, 
it  is  a  natural  extension  of  the  commutativity  property  in  ID  [1].  Note  that  in  such  case,  the  commutativity 
is  guaranteed  by  the  ID  Lie  group  structure  and  need  not  be  postulated.  Note  also  that  in  special  relativity, 
the  third  additional  postulate  is  isotropy.  This  leads  to  a  representation  in  term  of  rotations,  which  is  NOT 
commutative  in  D  >  1 .  There  is  therefore  no  hope  to  find  the  Lorentz  group  as  a  special  case  of  our  similarity 
transformations,  contrary  to  the  ID  case. 


With  this  additional  postulate,  the  form  of  the  matrix  (3.2)  can  be  completely  specified.  The  proof  is 
given  in  Appendix  A.  It  is: 

A  =  {l~-aU-jV)T  B  =  -{jU  +  dV)T 

(3.3)  D  ^  -c{jU  +  dV)T  E  =  {l-cdU-  bV)T 

G  ==  -hT  H  =  -kT 

where  F  is  a  function  of  U  and  V  obeying  the  transformation  rule: 

(3.4)  T{u\r^)  =  T(u\v')T{u,v){i  -  huu'  -  kvr), 

with  following  (3.7).  We  were  not  able  to  find  the  exact  expression  of  F  in  the  general  case, 

but  in  the  simple  cases  we  study  in  the  next  two  sections  (where  h  =  fe  =  0),F  =  lis  always  a  solution.  In 
general,  we  postulate  that  F  is  a  quadratic  form  in  U  and  F,  but  the  explicit  value  of  the  coefficient  is  not 
needed  in  the  derivation  of  the  composition  law. 
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The  coefficients  a,  6,  c,  d,  h,  k  and  j  are  parameters  depending  on  the  physical  model.  They  satisfy 


(3.5)  h  =  ck,  -  k ad  + jb  =  +  d^c 

They  are  real  if  the  scale  symmetry  is  continuous  [1]. 

Composing  two  transformations  from  reference  fields  R  =  (n,p)  to  R'  —  and  from  R'  = 

to  i?''  =  we  obtain: 

(3.6)  {UR\R>r,VjllRn)  =  (URlRf  ,VR\Rt)  {Ur>  iRlf  ,Vr> \Rrr) 
where  0  is  the  composition  law: 

-  aUU'  -  dVr  -  j{U^V  +  UV') 

1  -  hUU'  -  kVV'  , 

^  cjUU^  -  hVV‘  -  cd{WV  +  Ur) 

1  hUW  -  kVV' 

Like  in  ID,  the  composition  law  can  be  characterized  by  its  fixed  points,  which  depend  on  the  parameters 
6, /?,  j.  It  can  also  be  used  to  determine  the  scaling  exponents,  by  an  appropriate  choice  of  R,  R! 
and  R".  In  both  cases,  the  general  computation  is  rather  cumbersome  and  complicated  by  the  large  number 
of  degenerate  situations  (collapsing  or  vanishing  fixed  points).  We  treat  in  the  following  sections  two  special 
cases  where  the  computations  are  tractable.  We  emphasize  that  the  case  c  =  0  implies  no  influence  of  X 
over  Y  and  is  relevant  to  the  case  of  the  passive  scalar  dynamics  when  X  is  taken  as  the  passive  scalar 
log-coordinate  and  Y  as  the  velocity  field  log-coordinate. 

3.3.  The  choice  of  the  reference  field.  To  obtain  the  composition  law  for  the  scaling  exponents,  we 
must  make  precise  the  link  between  X,  Y  and  C/,  V.  This  depends  on  the  choice  of  <^o  ■A.s  emphasized 

in  [1],  a  convenient  choice  is  to  take  (po  =<  ln(<^^)  >  and  =<  ffi(V’^)  >.  Physically,  these  quantities  are 
the  most  probable  values  of  the  fields  (j)  and  The  log  coordinate  X  and  Y  then  characterize  deviations 
with  respect  to  these  most  probable  values,  and  are  a  measure  of  the  intermittency  of  the  fields.  In  such 
simple  case,  it  is  easy  to  check  [6]  that  the  scaling  exponents  for  0  and  -0  are  simply  connected  to  X  and  Y , 
via: 

dn^{n)  -  9„^(n)|„=o  = 

(3.8)  dpCip)  -  5pC(p)|p=o  =  ~- 

Now,  the  connection  between  U,V  and  ^(n),C(p)  can  be  made  more  clearly  by  exploiting  the  remark 
made  in  Section  III. A.  Indeed,  considering  the  transformation  from  X  =  0,y  =  0  to  =  (d„^(n)  — 
a„^(n)l„=o)T,yp  =  (dpCip)  -  dpC{p)\p=o)T,  we  get: 

U{n)  =  dn^{n)  -  a„^(n)|„=o  = 

dY 

(3.9)  V{p)=  dpC(p)  -  dpCip) |p=o  =  . 

This  shows  the  existence  of  a  bijection  between  {U,V)  and  (n,p).  It  unambiguously  defines  an  internal 
composition  law  (g)  on  (n,p),  through  a  transport  of  the  group  structure  0: 

(C7,  V){{n,p)§>{n\p*)  =  (U,V){n,p)<^{U,V){n\p'), 

(n,p)®(n',p')  =  (U-\V-^){{U,V){n,p) 

(3.10)  m,V){n’,p')). 


iu" 

< 

V" 
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From  this  property,  we  can  deduce  that  the  shape  of  (§)  is  the  same  as  (3.7),  with  different  parameters, 
characterizing  the  fixed  points  in  the  (n,p)  space^.  We  have  (n",p")  =  with: 

n  +  n'  -  ann'  -  Spp'  -  ^(n'p  +  np') 


(3.11) 


P"  = 


1  —  rjnn^  —  XPP* 

p^pf  -  rjnn^  -  ppp^  -  T6{n'p  +  np') 


1  —  7/nn'  —  XPP^ 

where  the  parameters  characterize  the  fixed  points  of  d)  and  satisfy: 

(3.12)  'n-TX,  +  +  =  + 


3.4.  Computation  of  the  scaling  exponents.  The  correspondence  established  through  eq.  (3.10) 
enables  to  classify  the  statistics  of  scale  invariant  random  process,  i.e.,  of  the  possible  shapes  for  ^(n),  C(p), 
on  universal  grounds.  Consider  eq.  (3.10)  with  (n,p)  =  (n',p')  —  (1,1).  This  provides  a  recursive  equation 
which  allows  the  computation  of  [U^V)  for  any  integer  m,  and  then,  by  continuation,  on  any  real  number, 
namely: 

E{m)  =  (1, 1)0(1, 1)0...0(1, 1)  =  (1, 1)'“^, 

{U,V){E{m))  =  (l/(l),F(l))0...0(l^(l),ni)) 

(3.13)  =(f/(l),F(l))[”*l. 

Here  the  notation  [m]  (resp.  [m]),  stands  for  the  iterate  via  0  (resp.  0):  we  have  introduced  the 
notation  E{m)  —  {E\{m)^E2{rn))  to  generalize  the  exponentiation  m  to  real,  defined  by  using  infinitesimal 
n  and  p.  It  is  now  technically  possible  to  calculate  the  values  of  the  C(^)  ^(p)'  set  of  equations 

(3.7,  3.10,  3.11,  3.13)  is  a  mathematically  well-posed  problem.  Indeed,  inverting  (3.8),  we  get: 


dgn) 

dn 


dn 


+  /7(l)[^i^^)], 


(3.14)  ®  =  ^U.o  +  l-CDft-'"". 

an  dn 

This  provides  a  symbolic  representation  of  the  possible  shape  for  the  scaling  exponents,  depending  on  the 
fixed  points  of  the  composition  law.  As  an  example,  we  treat  below  some  tractable  examples,  where  the 
fixed  points  are  simple.  Since  we  are  interested  in  the  case  where  no  coupling  exist  from  the  field  (j)  onto  ^ 
(passive  scalar  limit),  we  see  that  this  imposes  c  =  r  =  0.  These  values  are  adopted  from  now  on. 


4.  Passive  scalar  advected  by  a  Gaussian  field.  We  consider  here  the  case  where  u  is  a  Gaussian 
velocity  field,  with  2n-th  order  structure  functions  scaling  like: 

(4.1)  {{Suif^)  =  {{u{x  +  e)-  u{x)f)  OC 

Because  u  is  Gaussian,  it  follows  the  simple  “normal”  scaling  property: 


(4.2)  C(2n)  =  nC(2). 

We  are  interested  in  the  scaling  properties  of  a  passive  scalar,  advected  by  such  a  velocity  field.  Specifically, 
we  are  looking  for  the  scaling  exponents  ^(2n)  of  the  2n**  order  structure  function  of  temperature  increments 
defined  as: 

(4.3)  <  {e{x  -h  ^)  -  0{x)f^  >- 

^Physically,  these  fixed  points  correspond  to  the  limit  beyond  which  the  moments  of  the  distribution  in  (j)  or 'ip  are  divergent, 
see  [6]. 
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In  the  Kraichnan  model  [10]  of  the  eq.(l.l)  where  u  is  delta  correlated  in  time,  a  number  of  exact  results 
have  been  recently  established.  Kraichnan  [10]  and  Gawedzki  and  Kupianen  [11]  have  shown  the  second 
order  scaling  exponent  ^(2)  is  related  to  the  exponent  ^(2)  of  the  second  order  structure  function  of  the 
velocity  field  by: 

(4.4)  ^(2)=2-C(2). 

For  scaling  exponents  of  higher  order  structure  functions,  both  perturbative  and  non-perturbative  results 
are  available.  In  the  limit  where  C(2)  —>■  0,  it  was  shown  by  Bernard  et  al.  [12]  that  the  scaling  exponents 
follow  the  “anomalous”  law: 

(4.5)  ^(2n)  =  nm  -  C(2)  ^  0, 

where  D  is  the  space  dimension.  In  the  limit  where  C(2)  2,  it  was  shown  by  [13,  14]  that  the  anomalous 

correction  vanish,  so  that: 

(4.6)  ^(2n)=  71^(2),  C(2)^2. 

Chertkov  et  al,  [15]  used  an  expansion  in  the  inverse  1/D  of  the  space  dimensionality  to  obtain: 

(4.7)  ^(2n)  <(2)  -  ^^-C(2),  ^  0. 

which  is  consistent  with  the  result  (4.5)  obtained  by  Bernard  et  al,  [12]. 

On  the  other  hand,  Kraichnan  [10]  used  the  linear  ansatz  for  the  molecular-diffusion  term  and  found  a 
non-perturbative  result: 

(4.8)  C(2n)  =  ^V4nD^(2)  +  (D  -  mr  -  ^(^  "  ^(2)). 

Note  that  this  result  is  inconsistent  with  the  C(2)  expansion  of  Bernard  et  al  [12]:  using  (4.4),  one  finds  that 
in  the  limit  C(2)  -)■  0,  the  exponents  C(2n)  -  n^(2)  given  by  (4.8)  tend  to  a  finite  limit  which  contradicts 
(4.5). 

We  finally  mention  the  exact  non-perturbative  calculation  of  the  anomalous  correction  to  fourth  scaling 
exponent  p4  =  2^(2)  —  ^(4)  obtained  by  Benzi  et  al.  [16]  in  the  case  of  a  random  shell-model  for  passive 
scalar.  This  discrete  toy  model  mimics  most  of  the  properties  of  the  equation  (1.1)  [17].  In  such  model,  the 
properties  of  the  scaling  exponent  are  independent  of  ultraviolet  and  infrared  boundary  conditions,  except 
in  the  case  where  ^(2)  ^  0  if  the  diffusive  scale  is  kept  fixed.  In  this  last  case,  diffusive  effects  destroy  the 
inertial  range  scaling,  and  p4  0,  like  in  the  models  of  [12]  or  [15],  in  the  limit  C(2)  — ^  2.  Chertkov  et  al 
also  shows  p4  =  4(2  -  ((2))/D  [18]. 

4.1.  The  scaling  exponents.  Let  us  now  examine  those  predictions  within  the  framework  of  the  scale 
invariant  theory  obtained  in  Section  3.  In  the  stationary  situation  and  in  the  non-diffusive  limit,  the  equation 
(1.1)  is  invariant  under  the  transformation  Sh,h’{^)-  This  shows  that,  for  the  range  of  scale  in  which  the 
non-diffusive  limit  is  meaningful  (the  so-called  inertial  range),  the  structure  functions  obey  power  laws  (4.3), 
where  the  scaling  exponents  ^  (n)  may  be  a  priori  computed  by  using  scale  invariant  theory. 

In  the  Gaussian  case,  all  the  moments  are  convergent.  The  convergence  of  all  moments  implies  that 
the  nonzero  fixed  points  of  the  composition  law  for  p  must  be  —  oo  or  oo.  That  is  possible  only  if  r  =  7y  = 
X  =  p  =  0,  Since  the  scaling  exponent  of  the  Gaussian  field  is  linear,  that  is  Ci'^n)  =  nC(2),  V{n)  has  to 
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be  identically  zero  according  to  (3.8).  With  c  =  h  =  0  as  imposed  by  the  passive  scalar  requirement  and 
b  =  k  —  0  required  by  the  linear  scaling  exponent,  the  relevant  composition  law  becomes: 


(4.9) 

and 


(4.10) 


U  +  U'-  aUU' 
V  +  V 


n  +  n'  —  ann' 


Note  that  it  is  compatible  with  the  requirement  ^(p)  =  0.  The  solution  of  eqs.  (4.9)  and  (4.10)  can  be 
found  by  iteration,  as  explained  in  Section  III.D.  We  now  introduce  the  parameter  13  via 


(4.11)  ^^l-al7(l). 

The  solution  U{p)  of  the  equation  (4.9)  takes  two  different  forms  corresponding  to  the  two  following  cases: 

(1)  If  ^  ^  1  and  a  0,  then 


(4.12)  I7(n)  =  i  (l  -  (1  -  an)^“^) 

where  A  -  1  =  ln/9/ln(l  -  an(l)). 

Using  eq.(3.8)  and  integrating  over  n,  we  obtain: 

(4.13)  ^(2n)  =  2nho  “  ^  , 

where  /lo  is  a  constant  depending  on  5n^(0).  This  expression  is  a  slight  generalization  of  Kraichnan’s  formula. 
Indeed,  setting  /iq  =  0,  A  =  1/2,  o  =  -2{D  -  i{2))ID^{2)  and  a  =  -2D^{2)I[D  -  i{2)f,  we  find  exactly 
(4.8).  The  ’’log-normal”  formula  of  Chertkov  et  al  (4.7)  and  Bernard  et  al  (4.5),  with  ^(2n)  being  the  sum 
of  a  linear  and  a  parabolic  contribution  is  then  obtained  from  (4.13)  by  taking  the  limit  either  A  — >•  2,  or  the 
limit  as  follows. 

(2)  /?  =  1  and  a  =  0.  Then, 


(4.14)  ?(2n.)  n^(2)  —  2  n(n  1). 

Note  that  this  limit  is  not  achieved  when  the  parameters  are  “tuned”  to  Kraichnan’s  value  in  the  limit 
({2)  0,  ^(2)  2.  This  illustrates  in  another  way  the  discrepancy  arises  between  Kraichnan’s  result  and 

Bernard  et  aJ.’s  asymptotic  result. 

Finally,  we  get  two  additional  interesting  possibilities  when 

(3)  ^  —  1  and  a  /  0 


(4.15) 


17(1) 


in(l  —  an(l)) 
^(2n)  =  2nho  + 

U{1) 


ln(l  —  an), 


aln(l  —  an(l)) 


[2an  -  (1  -  2an)  ln(l  -  2an)]. 
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This  expression  tends  towards  the  log-normal  asymptotic  form  when  a  0. 
(4)  ^  1  and  a  =  0 


(4.16)  ^{2n)  =  2nho  +  ^  (l  -  . 

This  is  the  log-Poisson  law  [19,  20,  21],  which  again  tends  to  a  log-normal  in  the  limit  l3  1. 

4.2.  Discussion.  Our  result  obviously  generalizes  the  results  of  [10,  12,  15],  using  four  arbitrary  con¬ 
stants.  Because  we  use  symmetry  arguments,  we  cannot  obtain  any  quantitative  information  about  the 
value  of  the  anomalous  correction.  Only  “matching”  with  exact  computations  (starting  from  the  original 
equations)  can  provide  values  for  these  constants.  In  the  present  case,  we  have  four  arbitrary  parameters:  a, 
U{l)ya  and  /iq-  These  parameters  can  be  arbitrary  function  of  C{2)  which  need  to  be  determined.  The  exact 
result  (4.4)  enables  to  write  an  equation  linking  our  four  parameters  to  C(2).  Besides  that,  the  asymptotic 
expressions  of  Bernard  et  aJ.  and  Chertkov  et  al.  only  provide  asymptotic  expressions  for  ?7(1),  and  the 
condition  that  a  must  tend  to  zero  as  D  ->  oo  or  C(2)  0.  This  obviously  does  not  constrain  our  parameters 

enough,  so  we  shall  not  bother  to  try  to  exhibit  one  out  of  the  many  solutions  which  fits  all  the  requirement. 
It  would  be  interesting,  as  one  gets  more  and  more  precise  numerical  result  as  those  obtained  by  Frisch  et 
al.  [22],  to  try  to  derive  the  four  parameters  for  any  given  value  of  C(2),  and  see  if  they  are  enough  to  fit  the 
whole  collection  of  exponents  at  this  value  of  C(2).  This  would  provide  a  test  of  the  present  theory. 

5.  Passive  scalar  advected  by  turbulent  flow.  In  this  Section,  we  consider  the  case  where  the 
velocity  field  u  follows  the  Navier-Stokes  equations: 

VP 

(5.1)  dtVL  -h  (u  •  V)u  =  — —  +  uAv  -h  f. 

Here  p  is  fluid  density,  f  force,  P  pressure  and  v  viscosity.  Eq.  (1.1)  indicates  that  the  passive  scalar  is 
transported  by  the  turbulence.  Eq.  (5.1)  shows  that  the  dynamics  of  the  velocity  field  is  not  influenced  by 
the  passive  scalars. 

In  the  stationary  unforced  inviscid  limit  (i/,  k  ->  0),  the  set  of  equations  (1.1,5. 1)  is  invariant  under 
the  transformation  <S/i,/i/(A).  This  shows  that  for  the  range  of  scale  in  which  the  inviscid  unforced  limit  is 
meaningful  (the  so-called  inertial  range)  the  structure  functions  of  velocity  and  passive  scalar  obey  power 
laws  (1.2).  Dimensional  considerations  lead  to  the  prediction  ^(n)  =  (^{n)  ~  n/3-  However,  the  experiments 
and  numerical  simulations  do  not  support  the  predictions.  These  exponents  are  observed  to  deviate  from 
the  linear  laws  and  the  deviation  is  higher  for  the  temperature  than  for  the  velocity.  It  is  believed  that  the 
deviations  are  due  to  intermittency  of  velocity  and  scalar  dissipations. 

Many  models  have  been  proposed  to  reproduce  the  scaling  exponents  of  velocity  structures.  Among 
them,  She  and  Leveque’s  model  [19]  gives  an  excellent  prediction  of  experimental  results  without  any  free 
parameters.  Later,  this  model  is  shown  to  be  equivalent  to  assume  that  velocity  increments  follow  a  log- 
Poisson  statistics  [20,  23].  The  log-Poisson  statistics  converges  to  a  log-normal  statistics  under  a  suitable 
limit.  In  fact,  a  recent  careful  analysis  [24]  of  turbulent  data  using  wavelet  transform  indicated  that  the 
turbulence  is  so  close  to  log-normal  that  distinguishing  between  log-Poisson  and  log-normal  statistics  is 
impossible,  given  the  available  amount  of  turbulent  data.  Therefore,  we  shall  assume  in  the  following 
discussion  that  the  turbulent  velocity  is  log-Poisson,  and  then  take  the  suitable  limit  towards  log-normal 
statistics  to  examine  whether  the  passive  scalar  statistics  can  help  in  making  the  distinction. 
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As  an  extension  of  log-Poisson  statistics  of  velocity  fields,  Cao  and  Chen  [25]  assumed  a  bivariable  log- 
Poisson  distribution  for  joint  moments  of  velocity  dissipation  a  and  scalar  dissipation  This  assumption 
implies  that  passive  scalars  also  obey  log-Poisson  statistics,  which  is  in  remarkable  agreement  with  the  results 
of  numerical  simulations.  Ruiz,  Baudet  and  Ciliberto’s  recent  experiments  [26]  also  provide  strong  evidence 
to  support  log-Poisson  statistics  of  passive  scalars.  He  et  al.  [21]  propose  a  hierarchical  relation  for  the  joint 
moments  to  explain  why  passive  scalars  must  be  log-Poisson.  In  the  following,  we  explore  the  possibility  of 
log-Poisson  statistics  of  passive  scalars  within  the  formalism  developed  in  Section  3. 

5.1.  The  scaling  exponents.  A  remarkable  simplification  occurs  from  dynamical  considerations:  since 
velocity  dynamics  is  not  infiuenced  by  the  scalar  dynamics,  that  is,  the  composition  law  for  V  must  be 
independent  of  U,  one  must  choose  c  =  h  =  T  =  r]  —  0.  Moreover,  to  select  the  log-Poisson  statistics  for 
velocity  increments  [4]  requires  /c  =  0  and  j3  =  k  =  0.  Therefore,  the  relevant  composition  law  in  this  case  is: 


(5.2) 


and 


=  U' ’j- U  -  aUV  -  dV'V  -  jIU'V  +  UV'] 
V"  =  V'-hV-bV'V 


(5.3) 


n"  =  n'  -(-  n  —  an'n  —  Sp'p  —  7[n'p  +  np'] 
p"  =  p'  +  p 


A  rescaling  of  U  and  F  by  an  arbitrary  factor  t,  and  of  a,  6,  d  and  j  by  1/t  leaves  (5.2)  unchanged.  At 
the  same  time,  a  rescaling  of  n  and  p  by  1/t,  and  a  rescaling  of  a,  5,  7  by  a  factor  t  leaves  (5.3)  unchanged. 
Since  the  scaling  exponents  depends  only  on  the  derivative  of  U  and  V  with  respect  to  n  and  p,  we  may 
assume  without  loss  of  generality  that,  say,  p(l)  =  1. 

The  equation  (5.2)  has  one  stable  fixed  point 


and  one  unstable  fixed  point 

(5.5)  =  K.i 

We  solve  (5.2)  and  (5.3)  and  get  an  implicit  representation  of  the  solution  via. 


(5.6) 

where 


(5.7) 

and 


U[n{m)]  =  {U+-U-){\-pr)  +  y*(X-P^) 

V\pini)]  = 


A  =  l- 


U+-U-. 


Ih  =  l- 


K  • 
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(5.8) 

1  n(7n) 

1  vim) 

\  1  -  Oi  /  1  -  Oi  1  -  01 

=  m 

where 

Oi  =  1  -  an(l)  -  7, 

(5.9) 

02  =  ~  7n(l). 

Here,  n(l)  is  a  free  parameters.  To  avoid  any  discontinuity  while  ensuring  n(0)  =  0,  we  shall  only  assume 
that  n(0)  tends  to  when  ai  tends  to  0.  Note  here  that  the  log-normal  shape  for  the  velocity  increments 
is  obtained  by  letting  /32  ^  1,  in  which  case  the  function  U  simply  becomes: 


(5.10)  U{m)  =  {U^  -  U-)  (1  -  /3r)  +  rn. 

To  get  U{n)  as  a  closed  form,  one  needs  to  invert  the  function  n(m)  given  by  (5.8).  This  is  analytically 
possible  only  for  special  values  of  the  parameters,  otherwise  one  can  only  use  the  implicit  representation 
(5,6)  and  (5.8).  The  special  values  where  an  explicit  formulation  of  U{n)  is  possible  are: 

1)  when  02  =  0:  we  get: 

U{n)  =  {U+  -  U_)  (1  -  (1  +  n/a3)^'-i) 

+  T4{l-(l  +  n/a3)^="^), 

^{n)  =  nh^  +  . —  (l  -  (1  +  n/o3)'^‘) 

(5.11)  C(p)=Mc- j^(i-i3?), 

where  Aj  -  1  =  ln/?j/ln(ai),  for  i  =  1,2,  03  =  n(l)/(ai  -  1)  and  and  are  integration  constants.  The 
shape  for  C  is  the  usual  log-Poisson  shape.  The  shape  for  ^  is  a  superposition  of  log-levy  function.  In  the 
log-normal  limit  for  the  velocity  increments,  one  gets  a  different  shape: 

f  (n)  =  n/i£  +  (1  -  (1  +  n/03)^' ) 

-h  03  ((1  -h  n/03)  ln(l  -h  71/03)  -  n/03) , 

(5.12)  C(p)=p/ic  +  Y- 

So  the  log-normal  limit  is  not  really  simpler  than  the  general  solution. 

2)  when  oi  =  0,  we  get: 

U{n)  =  {U+  -  U-)  (1  -  +  K  (1  -  , 

V[p)  =  V.(l-pl), 

an)  =  nh^  -  (1  - 

(■-«'"'")• 

(5.13)  C(p)=pft(-rJ^(l -/??), 
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where  and  are  integration  constants.  In  this  case,  they  are  scaling  exponents  connected  with  the 
properties  of  the  rarest  events  (the  most  intermittent  structures),  whose  connexion  with  finite  size  effects  is 
discussed  in  [5].  The  physical  interpretation  of  the  parameters  appearing  in  (5.13)  is  obtained  by  taking  the 
limit  71  00  or  p  00,  when  /?i  and  132  are  less  than  one.  We  find  that 

^(n)  -n(l)(C/+-  ?7.-)/lnA  -n(l)K/ln^2 


and 


c(p)  ocp/i^  -K/ln/92. 


This  suggests  interpreting  —7i(l)({7_j.  —  UJ)l\nj3i  —  7i(l)K/ln/52  and  —V^l\np2  as  the  codimensions  of 
the  most  intermittent  structures  of  the  passive  scalar  and  of  the  velocity.  At  this  stage,  we  can  make  an 
interesting  remark:  when  n(l)  is  small  (this  is  a  case  of  weak  coupling  between  passive  scalar  and  velocity), 
the  last  term  in  (5.13c)  can  be  neglected.  In  this  case,  the  scaling  exponents  of  the  passive  scalar  have  a  log- 
Poisson  shape.  We  show  later  by  fitting  the  experimental  data  that  this  coupling  coefficient  is  indeed  small. 
This  therefore  explains  the  observation  in  [26].  Also,  we  note  that  the  scaling  exponents  (5.13)  reassemble  a 
log-Poisson  distribution  with  two  “atoms”,  whose  existence  was  conjectured  in  [23].  Finally,  it  is  interesting 
to  consider  the  limit  where  the  velocity  increments  are  log-normal.  In  that  case,  we  get: 


(5.14) 


^(n)  =  nh^  - 

Cip)  + 


n(l)(I/+-I/-) 
In  A 

2  ’ 


27z(1)  ’ 


In  that  case,  the  scaling  exponents  are  the  superposition  of  an  exponential  log-Poisson  shape,  and  of  a 
parabolic  log-normal  shape.  At  large  positive  n  (if  <  1),  the  parabolic  shape  dominates  the  scaling 
behavior,  and  one  does  not  observe  the  simple  linear  asymptotics  existing  in  the  log-Poisson  case.  This 
observation  could  be  used  as  a  tool  to  discriminate  between  log-Poisson  and  log-normal  velocity  increments, 
by  considering  the  statistics  of  passive  scalars.  This  is  an  indirect  check,  but  in  some  experimental  settings, 
this  might  actually  prove  easier  than  the  direct  check. 

In  all  other  cases,  we  may  only  obtain  an  implicit  representation  of  the  scaling  exponent  by  using  the 
chain  rule  d^fdm  =  d^/dn-  dn/dm. 


5.2.  Value  of  the  constants.  We  obtained  different  explicit  and  implicit  forms  of  scaling  exponents 
depending  on  several  constants.  In  the  most  general  (implicit)  case,  we  have  9  free  parameters.  In  the 
simplest  case  (number  2),  we  have  6  free  parameters.  We  need  the  same  number  of  constraints  to  compute 
these  parameters.  For  /3i  and  V*/ln/?27  we  can  use  the  values  obtained  by  [19],  which  provide  a  good 
fit  of  the  velocity  increments:  =  1/9,  /32  =  0.88  and  K/ln^2  =  2.  For  the  remaining  parameters, 

we  need  to  use  values  measured  for  the  passive  scalar.  Given  the  available  precision  of  the  measures,  it 
is  not  really  realistic  to  consider  more  than  n  =  5  or  72  =  6  scaling  exponents  for  computing  the  values 
of  the  constants.  It  is  then  only  possible  to  determine  the  constants  in  the  log-Poisson  case,  where  only 
three  remaining  parameters  need  to  be  determined.  This  simplification  enables  a  partial  determination 
of  the  parameters.  If  the  most  intermittent  structure  have  respectively  vortex  filaments  and  temperature 
sheets,  as  observed  in  numerical  simulations  [27,  28],  their  codimensions  are  approximatively  one  and  two. 
By  considering  that  the  most  intense  energy  dissipation  events  are  only  driven  by  the  mean  transfer  term 
(<  (5u)^  >)  and  using  the  K62  refined  similarity  hypothesis  [29],  She  and  Leveque  [19]  obtain  /i^  =  1/9,  and 
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Fig.  5-1.  Comparison  between  the  experimental  data  of  Ruiz  et  al.  (diamonds)  and  formula  (49)  with:  A  set  of  parameter 
derived  by  physical  arguments  copying  the  argumentation  of  She  and  Leveque  [19]  for  the  velocity  increments.  This  gives 
(/i^,/0i,n(l),C/+  -  U-,h^,V*,02)  =  (1/9,0.985,0.1,0.015,1/9,0.26,0.88)  (dash  line);  A  set  of  parameters  derived  by  a  best  fit 
to  the  data  {h^,0i,n(l),U+  -  =  (0.06,0.915,0.25,0.12,1/9,0.26,0.88)  (solid  line). 


then,  using  ((3)  =  1,  =  2/3.  A  same  kind  of  argument  also  leads  to  =  1/9.  Indeed,  if  we  suppose  that 

the  most  intense  temperature  dissipation  events  are  only  driven  by  <  >  and  <  Sv{S0)‘^  >,  we  obtain 

rv.  02  <  (Jv)3  >1/3  !  (Ji;(J0)2  where  0  is  some  constant  temperature.  Since  <  (<5?;)^  >~  i  and 

<  6v{60y‘  we  obtain  Then,  by  using  the  refined  Obukhov-Corrsin  hypothesis  ([30]),  we 

get  =  1/9  +  /iiv/2  + 1/3  =  1/9.  We  may  then  use  these  values  and  vary  n(l)  and  [ii  to  obtain  the  best  fit 
with  the  data.  The  result  is  given  in  Figure  1  for  n(l)  =  0.1  and  /3i  =0.86.  It  is  not  really  good,  especially 
at  large  n  where  the  formula  gives  a  much  larger  value  than  the  data.  This  is  because  the  value  =  1/9  is 
too  high.  A  three  parameter  fit  of  the  data  gives  =  0.06,  n(l)  =  0.25,  /?i  =  0.7,  resulting  in  a  codimension 
of  0.84  for  the  passive  scalar  rarest  events.  The  result  is  displayed  in  Figure  1. 

Certainly,  more  work  is  needed  to  understand  the  meaning  of  the  parameters.  A  recent  work  of  Dubrulle 
[5]  [31]  suggest  that  some  of  them  might  be  linked  with  finite  size  effects.  Work  is  in  progress  in  that 
direction. 

6.  Conclusion  and  discussion.  In  this  paper,  we  develop  a  formalism  leading  to  the  equations  fol¬ 
lowed  by  the  scaling  exponents  of  two-dimensional  scale  invariant  random  fields.  It  is  the  extension  of  [1]  on 
one-dimensional  random  fields.  We  do  not  make  any  assumption  on  (statistical)  relations  of  the  two  random 
fields.  For  example,  they  may  be  independent  of  each  other  or  anisotropic.  Thus,  the  two  random  fields  pos¬ 
sess  independent  scaling  exponents  and  they  may  have  different  limit  velocities.  This  situation  is  completely 
different  from  the  relativity  theory  of  motion,  where  the  limit  velocity  is  the  same  one  (light  velocity)  in 
each  direction.  We  begin  with  the  structures  of  transformation  group,  and  then  compute  the  transformation 
matrix.  The  composition  law,  derived  from  the  similarity  transformation,  indicates  the  equations  relating 
the  different  scaling  exponents.  Based  on  the  scaling  exponent  equations,  the  classification  of  all  exponents 
is  possible  by  the  computation  of  the  fixed  points  of  the  transformation.  This  is  performed  in  a  simple  case, 
the  passive  scalar  advected  by  i)  a  random  Gaussian  velocity  field;  and  ii)  a  turbulent  velocity  field. 
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In  the  Gaussian  case,  we  show  that  the  passive  scalar  increments  follow  a  log-levy  distribution  gener¬ 
alizing  the  solution  of  Kraichnan.  This  distribution  tends  in  an  appropriate  asymptotic  limit,  towards  a 
log-normal  distribution.  The  anomalous  correction  to  the  scaling  exponent  can  be  matched  with  the  Kraich¬ 
nan  solution  [10]  or  with  recent  perturbative  expansions  of  Chertkov  et  aJ.  [15]  and  Bernard  et  ai.  [12]  via 
an  appropriate  tuning  of  the  parameters,  which  cannot  be  predicted  by  our  simple  symmetry  arguments. 

In  the  turbulent  case,  we  show  that  when  the  velocity  increments  follow  log-Poisson  statistics,  the  passive 
scalar  increments  follow  various  different  statistics,  the  simplest  one  being  close  to  log-Poisson.  It  is  in  fact  a 
log-Poisson  distribution  with  two  “atoms”,  thereby  providing  an  explicit  example  of  a  random  system  with 
more  than  one  atom  whose  existence  was  conjectured  in  [23].  This  result  provides  an  explanation  of  the 
experimental  findings  of  [26]  about  the  temperature  increments. 

We  note  that  scale  invariance  may  play  a  role  in  many  other  physical  problems.  It  would  then  be 
interesting  to  consider  other  values  of  the  parameters  to  describe  other  physical  situations  such  as  scaling 
laws  in  magneto-hydrodynamics. 
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of  motion.  He  also  would  like  to  thank  SAP  for  its  hospitality.  B.  Dubrulle  acknowledges  useful  discussions 
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Appendix  A:  Derivation  of  Similarity  Transformations 

We  show  that  the  similarity  transformation  compatible  with  the  following  postulates  are  the  solution 
(3.3)  by  generalizing  the  method  [1]  to  the  two-dimensional  case: 

Posulate  1.  The  log-coordinates  are  invariant  under  global  translation. 

Posulate  2:  Among  all  imaginable  reference  fields,  there  exist  a  continuous  class  of  equivalent  reference 
fields  built  on  scale  invariant  fields. 

Posulate  3:  the  group  of  similarity  transformations  is  commutative. 

Postulate  2  implies  that  the  composition  of  two  similarity  transformations  with  respective  parameters 
(U,V)  and  should  be  a  similarity  transformation,  associated  to  a  third  parameter  noted 

Composition  of  two  transformations  in  (3.2)  then  gives 

A"  =  AA'+PD'-hPGW' 

B"  =  AB'  +  BE'  -h  TH'UV' 

r"u"  =  AV'u'  4-  BF'v'  +  rru 

D"  =  DA'  +  ED'aTG'VU' 

(6.1)  E"  =  DB'  aEE'  ATH'VV' 

r'v”  =  DT'u'  A  ET'v'  +  rrv 

G”U"  =  GA'U  +  HD'V  +  TG'U' 

H"V"  =  GB'U’VHE'V  +  TH'V' 

T"  =  GT'UU'  +  HV'VV'  -h  rr' 

We  now  use  the  Postulate  3  about  commutativity.  We  first  apply  the  transformation  {U' ^V')  and  then 
the  transformation  {U^V).  By  comparison  with  (6.1),  we  obtain  several  equations.  For  example,  equations 

(6.1.1)  and  (6.1.5)  give 

BD'-^B'D  =  UU'{T'G-TG'), 

^  ^  DB'-D'B  =  VV'{rH-TH'). 
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Factorizing  (6.2)  leads  to 
(6.3) 


D=cB,G  =  -hT,H  =  -kT, 


where  c,  h,  k  are  adjustable  parameters. 

Substituting  (6.3)  into  (6.1.2)  and  (6.1.4),  and  then  comparing  their  coefficients,  we  obtain 


h  =  ck. 


The  same  technique  applied  to  (6.1.3)  and  (6.1.6)  gives 

A  =  {l-^aU-jV)V 

(6.4)  B  =  -{jU  +  dV)T 

E  =  {l-cdU~hV)r 

where  a,b,d,j  are  adjustable  parameters.  (6.3)  and  (6.4)  are  solutions  of  equations  (6.1)  provided  the 
parameters  follow 

-k-\-ad-\-  jb  =  p  + 

This  completes  our  proof. 
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